i/i 


STUDY  OF^hT  CRITICAL  LAYER  IN  A  ROT^N(^UQUI[) 

“  PAYLQAD(U)  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT 

CENTER  ABERDEEN  PROVIN  .  R  SEDNEY  ET  AL.  AUG  84 
UNCLASSIFIED  ARBRL-TR-02582  S8I-AD-F300  478  F/G  1976 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS*  1963-A 


OTIC  FILE  copy  Agm,  AD- A 146  094 


AD  /V 


> 


» 


TECHNICAL  REPORT  ARBRL-TR-02582 

» 

A  STUDY  OF  THE  CRITICAL  LAYER  IN  A 
ROTATING  LIQUID  PAYLOAD 

i 


Raymond  Sedney 
Nathan  Gerber 


August  1984 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER 

BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


» 


Approved  for  public  release;  distribution  unlimited. 


I 


84  09  27  025 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.  S.  Department  of  Commerce,  Springfield,  Virginia 


22161. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position,  unless  so  designated  by  other 
authorized  documents. 


The  use  of  trade  names  or  manufacturers'  names  in  this  report 
does  not  constitute  indorsement  of  any  commercial  product. 


SECURITY  CLASS! FI  CATION  OF  THIS  PAGE  (Whon  Dt 


i  tormd) 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 


TECHNICAL  REPORT  ARBRL-TR-02582 


4.  TITLE  (and  Subtltlm) 

A  STUDY  OF  THE  CRITICAL  LAYER  IN  A 
ROTATING  LIQUID  PAYLOAD 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  4  PERIOD  COVEREO 


Final 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORS  A) 

Raymond  Sedney 
Nathan  Gerber 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Ballistic  Research  Laboratory 
ATTN:  URXBR-LF3 

Aberdeen  Proving  Ground,  Maryland  21005-5066 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Ballistic  Research  Laboratory 
ATTN:  DRXBR-OD-ST 

Aberdeen  Proving  Ground,  Maryland  21005-5066 


MONITORING  AGENCY  NAME  4  ADDRESSES  dlltaront  from  Controlling  Ollier) 


14.  DISTRIBUTION  STATEMENT  (of  thl.  RoporQ 


8.  CONTRACT  OR  GRANT  NUMBERf..) 


tO.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 


RDT&E  1L161102AH43 


12.  REPORT  DATE 

I  August  1984 

.  13.  NUMBER  OF  PAGES 


15*.  DECL  ASSl  F|  CATION  /  DOWNGRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  tho  obmtrmct  mntormd  in  Block  20 ,  If  diffmront  from  Report) 


19.  KEY  WOROS  (Contlnuo  on  tmvmtmm  oldo  it  nmcoomory  mtd  Identity  by  block  number) 


Critical  Layer 

Eigenvalue  Problem  -  Rotating  Flow 
Liquid-Filled  Projectile 
Rotating  Fluid 


Spin-U 


20.  AMT* ACT  fCanitmto  on  rov 


I  Idonilfy  by  block  rmmbor) 


(bja) 


The  unsteady  motion  of  a  fluid  which  fills  a  spinning  cylinder  is  con¬ 
sidered.  Spin-up  of  the  fluid  is  the  basic  flow  which  is  perturbed.  Non- 
ax  i  symmetric,  viscous  perturbations  are  used  to  study  the  wave  motion  and 
critical  layer.  The  frequencies  and  decay  rates  are  determined  by  the  eigen¬ 
values  of  the  system  of  perturbation  equations.  For  large  time  the  fluid 

(continued) 


corn  or  or  t  nov  si  is  obsolete 


_ UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Sim  Oat*  Entrrrd) 


_ UNCLASSIFIED _ 

SSCUWITY  CLAgjlFICATtOM  OP  THU  PAOKflmm  DaU> 


approaches  solid  body  rotation;  for  this  state  there  is  no  critical  layer  and 
the  eigenvalue  problem  is  considerably  simpler.  The  critical  layer  always 
exists  for  small  time;  it  ceases  to  exist  at  a  time  which  depends  on  the 
parameters  of  the  basic  flow  and  the  wave  motion.  Time  histories  of  the  eigen' 
values  and  of  the  critical  layer  are  given  for  two  cases  and  two  radial  modes. 
The  effects  of  the  critical  layer  on  the  eigenfunctions  and  the  phase  of  the 
velocity  are  presented.  Comparison  with  experiment  is  discussed. 


UNCLASSIFIED 


SCCURITV  CLASSIFICATION  OP  THIS  PAOK(*fta)i  Data  Mnf 


TABLE  OF  CONTENTS 


LIST  OF  ILLUSTRATIONS . 

INTRODUCTION . 

THE  SPIN-UP  BASIC  FLOW.... 

THE  PERTURBED  FLOW . 

THE  EIGENVALUE  PROBLEM.... 

A.  The  Mathematical  Model 

B.  The  Numerical  Method.. 

THE  CRITICAL  LAYER . 

RESULTS . 

DISCUSSION . 

ACKNOWLEDGMENT . 

REFERENCES . 

DISTRIBUTION  LIST . 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1  V  vs  r  for  Cases  1  and  2  at  Three  Times;  k  -  0.5 .  23 

2  The  Magnitude  of  the  Characteristic  Determinant,  |z|, 

vs  Cr  for  Three  Values  of  Cj  for  Case  1,  m  =  1, 

k  *  3,  t  *  155.  Cr2  *  .072614,  CI2  =  .034266 .  23 

3  The  Eigenfunction  Real  (w)  vs  r  for  Re  =  5  x  105, 

A  *  2.679,  Mode  (5,1,1),  ts  *  3788.  Note  Scale 

Changes . 24 

4  Time  Histories  for  Cj  and  C2  for  Case  1,  Modes 

(3,1,1,)  and  (3,2,1),  t$  =  460,  <  *  0.5 .  24 

5  Critical  Level,  r_,  vs  t  for  Cases  1  and  2, 

n  =  1  and  2,  ie  =  0.5 .  25 

6  CR1  vs  t  for  Case  2,  Mode  (5,1,1)  and  for  Two 

Variations  on  Case  2 .  25 

7  The  Phase  Angle  of  u'  for  Case  2,  Mode  (5,1,1) 

at  Three  Times .  26 

8  In  the  (r,e)  Plane,  the  Approximate  Have  Surface  u'  « 

Constant  for  Fixed  z'  and  t  for  Case  2,  Mode  (5,1,1) 
at  t  =  6,000.  The  Circular  Arc  Portion  Is  at  the 
Critical  Level  r  =  rc  =  .20;  the  Change  In  e  Near 

r  -  1  Occurs  In  the  Boundary  Layer .  26 

9  The  Phase  Angle  of  u'  for  Case  2,  Mode  (5,2,1)  at 

t  =  6,000.  The  Critical  Level  Is  at  rc  *  .303 .  27 

10  The  Phase  Angle  of  u1  for  Case  1,  Modes  (3,1,1)  and 

(3,2,1)  at  t  *  160 .  27 

11  CR  vs  t:  Experimental  Points  and  Calculated  Values; 

rc  vs  t,  k  -  0.443 .  28 


I.  INTRODUCTION 


The  unsteady  motion  of  a  fluid  which  fills  a  spinning  right-circular 
cylinder  is  considered  in  this  paper.  The  spin  is  imparted  impulsively  to  the 
cylinder  and  the  spin-up  of  the  fluid  is  perturbed  to  stutty  the  wave  motion  in 
the  rotating  fluid.  This  is  called  the  spin-up  eigenvalue  problem;  abbrevi¬ 
ated  e.v.  problem.  Because  of  boundary  layers  and  the  critical  layer  that 
exist  in  the  flow,  the  problem  is  formulated  with  viscous  perturbations.  The 
physical  significance  of  the  critical  layer  is  discussed  as  well  as  how  it 
affects  the  solution  of  the  numerical  problem.  The  term  "critical  layer"  is 
used  here  because  it  is  analogous  to  the  critical  layer  that  occurs  in  the 
Orr-Sommerfeld  equation  governing  stability  of  a  laminar  shear  layer;  although 
the  stability  of  the  rotating  flow  is  not  a  question  here,  the  formulation  of 
the  e.v.  problem  is  mathematically  the  same  as  for  the  stability  problem.  For 
the  Orr-Sommerfeld  equation,  a  critical  layer  exists  if  the  disturbance  phase 
velocity  is  equal  to  the  basic  flow  shear  velocity.  For  the  rotating  fluid, 
the  condition  for  a  critical  layer  is  that  the  wave  frequency  be  an  integral 
multiple  of  the  local  angular  frequency  of  the  basic  circumferential  flow. 

The  application  of  this  work  is  to  the  stutjy  of  the  flight  of  liquid- 
filled  projectiles;  these  have  a  proclivity  for  unusual  flight  behavior,  often 
being  unstable  even  though  the  same  projectile  with  a  solid  payload  is 
stable.  A  knowledge  of  the  wave  motion  in  the  rotating  fluid  is  fundamental 
to  an  understanding  of  the  effect  of  the  liquid  on  the  projectile  motion. 

The  frequencies  and  decay  rates  of  the  waves  are  determined  by  the  com¬ 
plex  eigenvalues  of  the  system  of  perturbation  equations.  For  large  time  the 
fluid  approaches  solid  body  rotation;  for  this  state  there  Is  no  critical 
layer  and  the  eigenvalue  problem  Is  considerably  simpler.  The  critical  layer 
always  exists  for  small  time;  it  ceases  to  exist  at  a  time  which  depends  on 
the  parameters  of  the  basic  flow  and  the  wave  motion. 

The  physics  of  the  basic  flow,  spin-up  from  rest,  was  presented  by 
Wedemeyer.1  The  flow  Is  determined  by  Reynolds  number  »  Re  »  a  a2/v  and 
aspect  ratio  *  A  *  c/a  where  a  is  the  spin  (rad/  sec),  a  and  c  are  the  radius 
and  half-height  of  the  cylinder  and  v  is  the  kinematic  viscosity  of  the 
fluid.  The  model  determines  the  core  flow,  not  that  in  the  endwall  boundary 
layers.  The  flow  can  also  be  determined  in  a  more  complete  way  by  solving  the 
Navier-Stokes  equations  by  finite  difference  methods.  For  the  e.v.  problem  it 
would  be  Impractical  to  use  the  finite  difference  solution  for  the  basic  flow. 

A  discussion  of  previous  attempts  to  solve  the  e.v.  problem  is  given  in 
Reference  2  together  with  more  details  of  the  work  presented  here.  Experi¬ 
mental  data  such  as  velocity  fields,  pressures,  or  gyroscopic  motion  of  the 


2.  E.  H.  Wedemeyer,  ” The  Unsteady  Flow  Within  a  Spinning  Cylinder,"  Journal 
of  Fluid  Meohanioe.  Vol.  20,  Part  3 ,  1964,  pp.  383-399 .  (See  also  BRL 
Report  Vo.  1225,  October  1963 ,  AD  431846.) 

2.  R.  Sedney  and  N.  Gerber,  "Oeoillatione  of  a  Liquid  in  a  Rotating 
Cylinder:  Part  II.  Spin-Up,"  US  Army  Ballietio  Reeearoh  Laboratory, 
Aberdeen  Proving  Ground,  MD,  BRL  Teohnioal  Report  ARBRL-TR-02489 ,  May 
1983.  (AD  A129094) 
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container  exist  for  the  range  of  parameters  1  <  Re  <  107  and  1  <  A  <  5.  The 

spin-up  basic  flow  and  the  e.v.  analyses  are  restricted  to  the  higher  Re  by 
virtue  of  some  asymptotic  approximations.  The  lower  bound  on  Re  for  applica¬ 
bility  of  the  e.v.  results  Is  not  known.  In  general.  For  axisymmetric 
perturbations  It  appears  to  be  0(102);  see  Reference  3.  At  present  the 
analysis  of  the  rotating  fluid  problem,  per  se.  Is  also  limited  because  the 
boundary  layers  on  the  cylinder  endwalls,  i.e.,  the  Ekman  layers,  and  the 
Stewartson  layer  on  the  sidewall  are  not  Included.  For  application  to  the 
projectile  problem  the  angular  motion  must  be  restricted  to  small  angles 
because  the  theory  Is  linearized. 

Whether  or  not  the  effects  of  the  time  dependent,  spin-up  basic  flow  are 
Important  In  a  projectile  flight,  rather  than  solid  body  rotation  of  the 
fluid,  can  be  estimated  by  comparing  the  characteristic  time  for  spin-up,  €$, 

with  the  time  of  flight  of  the  projectile.  If  the  former  Is  small  compared  to 
the  latter,  spin-up  effects  can  be  neglected.  For  laminar  Ekman  layers, 

C$  *  (2c/a)  Re*/2/n  *  2/E*/2  n  (sec) 

where  E  =  v/fl  c2  Is  the  Ekman  number,  often  used  In  rotating  fluid  problems 
rather  than  Re.  This  estimate  Is  derivable  from  linear  spin-up  theory4*5  or 
the  Wedemeyer  model.  We  use  the  non-dimensional  characteristic  spin-up  time 
ts  *  a€$.  For  Re  >  105,  approximately,  the  Ekman  layers  may  be  turbulent,  in 

which  case  the  characteristic  spin-up  time  can  be  estimated  by 

Est  -  (28.6  c/a)  Re1/5/n  (sec) 


which  can  be  obtained  from  the  Wedemeyer  solution  without  diffusion  for  turbu¬ 
lent  Ekman  layers.1  These  times  do  not  measure  how  close  the  flow  is  to  solid 
body  rotation.  A  rule  of  thumb  often  used,  but  not  always  accurate.  Is  that 
solid  body  rotation  Is  reached  at  about  4E$  after  an  Impulsive  angular 

velocity  Is  applied  to  the  cylinder. 

To  appreciate  these  time  scales  for  projectile  applications,  consider  two 
cases  which  will  be  used  to  present  Illustrative  numerical  results: 


3,  R.  Sidney,  It,  Gerber,  and  J,  M,  Bartoe ,  "Oeoillatione  of  a  Liquid  in  a 
Rotating  Cylinder "  AIAA  20th  Aerospace  Sciences  Meeting,  January  11-14 , 

1982,  Orlando,  Florida,  AIAA  Paper  82-0296,  (See  also  ARBRL-TR-02488 ,  May 

1983,  AD  A129088.) 

4,  H,  P,  Greenepan,  The  Theory  of  Rotating  Fluids .  Cambridge  University 
Preee ,  London  and  Reu  York,  1988. 

5,  S,  R.  Benton  and  A,  Clark,  Jr,,  "Spin-Up,"  article  in  Annual  Review  of 
Fluid  Mechanics,  Vol,  8,  Annual  Reviews,  Inc,,  Palo  Alto,  CA,  1974, 
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Re 

c/a 

ft  (rad/sec) 

fs  (sec) 

Case  1: 

4,974 

3.30 

8937 

0.052 

Case  2: 

1.99x10® 

5.20 

754 

3.510 

These  parameters,  courtesy  of  Dr.  W.  P.  D'Amico,  are  appropriate,  for  Case  1, 
to  a  test  of  small  caliber  projectiles  in  a  ballistic  range,  and  for  Case  2  to 
an  artillery  projectile.  At  the  spin-up  times  the  projectiles  would  be  48  m 
and  1,230  m  from  the  gun  for  Cases  1  and  2,  respectively;  in  both  cases  obser¬ 
vations  on  projectile  motion  could  be  made  at  these  distances.  Partial 
validation  of  the  theory  has  been  provided  by  experiments6  and  numerical 
simulations. 3 

Some  of  the  first  part  of  this  report  repeats  what  is  in  Reference  2.  It 
is  included  here  for  convenience. 


II.  THE  SPIN-UP  BASIC  FLOW 

Consider  the  axisymmetric,  time  dependent  motion  of  a  fluid  which  fills  a 
cylinder,  initially  at  rest,  which  is  impulsively  brought  to  a  constant 

angular  velocity  ft  with  respect  to  its  axis.  In  practice,  an  impulsive  start 
is  impossible;  the  conditions  for  approximating  it  and  the  degree  of 

approximation  in  some  experimental  apparatus  and  in  a  gun  tube  are  discussed 
in  Reference  7.  For  an  impulsive  start,  the  validity  of  the  Wedemeyer  model 
can  be  discussed  in  terms  of  three  time  scales:  the  time  for  one  revolution 

of  the  cylinder,  2rrft-1,  t$  or  and  the  time  for  vorticity  to  diffuse 

radially.  Re  ft-1.  The  model  requires  2irft_1  «  «  Re  ft-1. 

It  is  known  that  the  Ekman  layers  form  and  become  essentially  steady  in 

time  2nft_1.  Although  Wedemeyer1  showed  the  crucial  importance  of  the  Ekman 
layers  to  the  spin-up  process,  his  model  did  not  require  a  solution  for  the 
flow  in  these  layers;  exclusion  of  this  solution  has  important  consequences 
for  the  e.v.  problem.  The  basic  mechanism  for  spin-up  starts  with  the  suction 
exerted  by  the  Ekman  layers  which  draws  external  fluid  into  them  where 

rotation  is  induced.  With  no  pressure  gradient  acting,  the  fluid  spirals  out 
to  larger  radii  where  the  layers  eject  fluid.  At  those  radii  the  flow  outside 
the  Ekman  layers,  called  the  core  flow,  is  now  rotating.  This  mechanism  is 


6.  S.  Stergiopouloa ,  "An  Experimental  Study  of  Inertial  Waves  in  a  Fluid 
Contained  in  a  Rotating  Cylindrical  Cavity  During  Spin-Up  from  Rest," 
Ph.D.  Thesis,  York  University ,  Toronto ,  Ontario ,  February  1982. 

7.  R.  Sedney  and  N.  Gerber ,  "Viscous  Effects  in  the  Wedemeyer  Model  of  Spin- 
Up  From  Rest,"  US  Army  Ballistic  Research  Laboratory ,  Aberdeen  Proving 
Ground ,  MD.,  BRL  Technical  Report  ARBRL-TR-02493 ,  June  1983 .  (AD  A129506) 


much  more  efficient  than  diffusion  which  plays  a  small  role  in  spin-up. 
Wedemeyer1  showed  in  his  remarkable  paper  that  the  flow  must  be  divided  into 
two  regions:  the  Ekman  layers  and  the  core  flow.  Actually  an  additional 
boundary  layer,  a  Stewartson  layer.  Is  required  at  the  side  wall. 

In  the  following,  lengths,  velocities,  pressure,  and  time  are  made 

2  2  -1 

nondimensional  by  a,  an,  pn  a  ,  and  n  ,  respectively,  where  p  is  the  liquid 
density.  In  the  inertial  frame  cylindrical  coordinates,  r,  e,  z  are  used, 
with  the  origin  of  z  at  the  center  of  the  cylinder,  and  velocities  are  U,  V 
W,  respectively.  Dimensionless  time  is  t.  Derivatives  are  indicated  by 
subscripts. 

Starting  from  the  Navier-Stokes  equations  for  axisymmetric  flow  Wedemeyer 
used  order  of  magnitude  arguments  to  simplify  them  in  the  core  flow.  The 
three  momentum  equations  reduce  to 


Vt  +  U  <Vf  +  V/r)  =  Re'1  [Vpr  +  (V/r)r]  (2.1) 

and 

Uz  =  Vz  =  Pz  =  0  (2.2) 

For  Re  +  •  he  proposed  neglecting  the  diffusion  terms  in  (2.1)  so  that 

vwt  ♦  uw  <vwr  ♦  V>  *  0  <2-3> 


where  the  subscript  w  indicates  this  approximation.  Wedemeyer  used  (2.3) 
rather  than  (2.1)  when  he  applied  his  model. 

To  solve  (2.1),  (2.3)  a  relationship  between  U  and  V  is  necessary. 
Wedemeyer  used  the  facts  that  Ekman  layers  are  steady  after  one  revolution  and 
that  the  radial  mass  flux  in  the  core  flow  must  be  balanced  by  that  in  the 
Ekman  layers  to  obtain  some  conditions  on  the  U,  V  relationship.  At  this 
point  he  was  forced  to  take  a  phenomenological  approach.  The  relationship  is 
known  at  t  -*■  0  and  t  •  and  he  proposed  a  linear  interpolation  between  them 
to  obtain  an  approximate  relationship  for  any  t.  He  tested  this  idea  in  some 
other  problems  where  the  solution  was  known  and  decided  it  was  satisfactory. 
Some  confusion  has  appeared  in  later  literature  because  this  step  was 
misinterpreted;  this  matter  is  discussed  in  Reference  7.  The  result  is 

U  =  k^V-r)  k*  =  <(a/c)Re"1/2  =  2</ts  (2.4) 


for  laminar  Ekman  layers.  Wedemeyer  proposed  <  =  0.443  but  Greenspan4 
suggested  <  =  0.5;  the  latter  often  gives  results  in  better  agreement  with 
numerical  solutions  to  the  Navier-Stokes  equations.  Other  relationships  have 
been  proposed,  but  they  will  not  be  discussed  here;  see  Reference  7.  For 
turbulent  Ekman  layers 
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where  tst  is  the  nondimensional ,  turbulent  spin-up  time.  The  core  flow  is 

assumed  to  be  laminar  so  that  turbulent  stresses  are  not  introduced  in  the 
right-hand  side  of  (2.1). 

Using  (2.4),  (2.3)  can  be  solved  explicitly: 


Vw  =  (r  e  -  l/r)/(e  1  -  1)  for  r>e  1  (2.6) 


'M 

and  Vw  *  0  otherwise.  Therefore,  r  =  rf  =  e  separates  rotating  and 

nonrotating  fluid  where  there  is  a  discontinuity  in  shear.  W  is  obtained  from 
the  continuity  equation:  W  =  -( z/r)  (rU)r.  At  r  =  1,  W  *  0;  thus,  the 
Stewartson  layer  should  be  included  at  r  =  1,  but  this  has  yet  to  be  done. 
Also  W  *  0  at  the  end  walls  z  =  ±  c/a. 

Wedemeyer  pointed  out  that  the  corner  in  the  solution  (2.6)  would  be 
smoothed  out  if  the  diffusion  terms  were  retained,  as  in  (2.1).  But  this  is  a 
nonlinear  second  order  equation  and  must  be  integrated  by  finite  difference 
methods.  This  can  be  done  with  some  standard  techniques  for  diffusion 
equations.  Special  treatment  is  needed  near  the  point  r  =  1,  t  =  0  because, 
for  an  impulsive  start,  a  discontinuity  in  the  boundary  conditions  exists  at 
that  point.  In  our  work  a  local,  analytic  solution  was  derived  to  resolve  the 
discontinuity;  see  Reference  8.  In  most  of  our  e.v.  calculations  we  used  the 
V  from  the  numerical  solution  of  (2.1)  with  either  (2.4)  or  (2.5),  but  other 
options  are  also  used.  For  laminar  Ekman  layers,  the  spin-up  velocity  profile 


V  =  f  (r,  k^t,  kfc  Re); 


for  the  turbulent  case  k^  is  replaced  by  kt. 

Some  examples  of  the  solutions  of  (2.1)  will  be  given  to  illustrate  the  V 
profiles  which  must  be  perturbed  in  the  e.v.  problem.  The  parameters  of  Cases 
1  and  2  are  used  in  V(r,t)  presented  in  Figure  1.  For  Case  1,  the  relation 


8 .  R.  Sedney  and  N.  Gerber ,  "Treatment  of  the  Discontinuity  in  the  Spin-Up 
Problem  with  Impulsive  Start ,”  US  Amy  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  HD. ,  BRL  Technical  Report  ARBRL-TR-02S20 , 
September  1983 .  (AD  A133682) 


for  laminar  Ekman  layers,  (2.4),  Is  used  and  for  Case  2,  (2.5)  Is  used.  For 
each  case  the  results  for  three  times  are  shown.  As  t  -*■  +  0,  V  tends  to  a 
discontinuous  function;  for  small  t,  an  asymptotic  solution  to  (2.1)  was 
derived  in  Reference  8  using  erf  functions.  The  e.v.  analysis  uses  a  quasi¬ 
steady  assumption  which  Is  violated  as  t  ♦  +  0.  A  practical  limit  on  how 
small  t  should  be  In  the  analysis  Is  set  by  the  smallest  t  at  which  data  can 
be  obtained. 


III.  THE  PERTURBED  FLOW 

The  procedure  for  obtaining  the  equations  governing  the  perturbed  flow  is 
a  standard  one  and  will  only  be  outlined.  The  velocity  components  and  pres¬ 
sure,  expressed  as  the  sum  of  the  spin-up  basic  flow  and  the  perturbation, 
e.g.,  U*  (r,z,t)  +  u'(r,  0,z,t),  are  substituted  into  the  Navier-Stokes 
equations  for  3-D,  unsteady  flow.  Initially  the  basic  flow,  e.g.,  U*,  is  a 
solution  to  the  Navier-Stokes  equations  for  axi symmetric  flow,  the  zeroth 
order  terms.  The  first  order  terms,  linear  in  the  perturbations,  are  retained 
and  the  second  order  terms  are  neglected.  The  coefficients  in  the  linear 
perturbation  equations  are  the  basic  flow  variables  and  their  derivatives. 
The  basic  flow  is  now  approximated  by  the  results  from  the  Wedemeyer  model. 
All  coefficients  containing  U*  and  W*  and  their  derivatives  with  respect  to  r 
and  z  are  0(1/0  -  except  one  term  discussed  below;  all  those  containing  V* 
are  0(1).  Replacing  V*  by  V,  the  resulting  perturbation  equations  contain 
only  V,  not  U  and  W;  they  are 


u'  +(V/r)u*  -2Vv'/r  =  -p‘  +Re”*(v2u'-u7r2-2v'  /r2) 
u  9  r  o 


v't+[Vr+(V/r)]u'+(V/r)v,0  =  -p'  g/r+Re”*^'  -  v7r2+2u'0/r2  (3.1) 


w't+(V/r)w'0  =  -p^+Re’^w' 


(ru')r+v'0+  rw'z  =  0. 


These  equations  govern  viscous  perturbations  of  the  core  flow.  If  a 

-1/2 

formal,  rational  expansion  with  Re  as  the  parameter  had  been  used  the 

viscous  or  Re”*  terms  would  not  appear.  These  govern  the  inviscid  per¬ 

turbations  and  are  adequate  except  in  the  perturbation  boundary  layer  and  the 
critical  layer.  The  viscous  terms  could  be  included  by  local  analyses  in 
these  two  regions.  We  have  Included  them  in  a  global  sense  by  retaining 

the  Re"*  terms.  In  the  derivation  of  (3.1),  one  term  in  the  z-momentum 
equations  contains  W*r.  If  this  is  approximated  by  Wp  that  term  is  not  small 
near  r  *  1  because  the  Wedemeyer  model  requires  a  Stewartson  layer  there,  as 
discussed  earlier.  Therefore,  strictly,  (3.1)  are  valid  outside  the 

Stewartson  layer;  this  restriction  does  not  appear  to  be  a  serious  one. 
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Finally,  employing  the  quasi-steady  assumption  in  (3.1)  is  considered.  The 
time  scale  for  the  perturbation  is  0(1);  that  for  V  is  0(ts)  where  t$  »  1, 

except  for  t  +  0  as  discussed  earlier.  With  that  exception  the  variation 
of  V  with  t  in  (3.1)  is  slow  and  t  can  be  regarded  as  a  parameter  in  V(r,t). 


IV.  THE  EIGENVALUE  PROBEM 
A.  The  Mathematical  Model. 

It  is  convenient  to  introduce  the  new  coordinate  z'  =  z  +  A.  It  is 
assumed  that  the  perturbation  can  be  represented  as  a  superposition  of 
modes.  With  the  quasi-steady  assumption  such  a  separation  of  variables  is 
possible.  In  complex  notation 

u'  =  Real  ju(r)  cos  Kz‘  exp  [i(Ct-me)]|  (4.1) 

with  similar  expressions  for  v'  and  p';  w'  has  the  same  form  except  for  a 
sin  Kz'  factor.  Here  K  =  kir/2A  with  k  =  1,2,...  and  m  =  0,  ±  l,...the  axial 
and  azimuthal  wave  numbers,  respectively;  m  =  1  is  the  value  relevant  to  the 
projectile  problem.  The  complex  quantities  u,v,w,p  are  solutions  of  the 
system  of  ordinary  differential  equations  obtained  by  substituting  the  modal 
forms  into  (3.1).  The  nondimensional  complex  constant 

C  =  +  i  Cj 

is  the  eigenvalue  of  the  system.  The  dimensional  wave  frequency  is  CRn  and 
the  decay  rate  is  Cjfl. 

For  the  free  oscillation  problem  the  boundary  conditions  at  the  end  walls 
z’  =  0,  2A  are  u*  ■  v'  *  w*  ®  0  if  the  complete  flow  is  being  perturbed. 
Since  we  are  perturbing  the  core  flow,  not  the  Ekman  layers,  the  boundary 
conditions  are  not  the  same.  The  modal  forms  give  w*  =  0  but  u'  and  v*  do  not 
satisfy  the  no-slip  condition.  Appropriate  inner  expansions  in  the  Ekman 
layers  are  required  to  correct  this.  Here  the  modal  form  is  used  without  end- 
wall  correction.  The  modal  form  is  also  used  for  the  solid  rotation  case  for 
which  there  are  no  Ekman  layers.  However,  endwall  corrections  are  still 
necessary;  the  perturbation  boundary  layer  must  be  introduced.  Not  including 
this  correction  gives  values  of  Cj  which  are  incorrect  by  a  factor  of  2  but 

the  effect  on  CR  is  about  1%.  An  ad  hoc  method  of  including  an  endwall 

correction  for  the  spin-up  case  is  to  use  the  solid  rotation  correction  for 
finite  time.  This  ad  hoc  correction  is  advantageous  for  large  t. 

From  (3.1)  a  sixth  order  system  of  ordinary  differential  equations  Is 
obtained: 


[Re”1(  Aj-r"2)-i  M]u+(2/r)(V+1  m  r"1Re"1)v-p(.  «  0 


[Re-1(A1-r"2)-i  M]v-( 3V/3r+V/r+2i  m  Re'^/r2^  +  1  m  p/r  *  0  (4.2) 


[Re^^-i  M]  w  +  K  p  =  0 


(r  u)p  -imv  +  Krw  =  0. 


where 

Alf  5  frr  +  fr/r  ’  £<n,2/r2)  + 


and 

M(r)  =  C  -  mV/r  (4.3) 


The  no-slip  boundary  conditions  at  r  *  1  require 


u*v*w  =  0  at  r  *  1.  (4.4) 


The  boundary  conditions  at  r  =  0  depend  on  m;  they  are  derived  using 
continuity  and  single-valuedness.  For  m  ■  1 


u-iv  =  w»p  =  0  at  r  »  0.  (4.5) 


The  system  (4.2),  (4.4),  and  (4.5)  is  not  self-adjoint.  Although  there 
is  no  proof,  we  assume  the  e.v.  form  a  denumerable,  discrete  spectrum.  With 
index  n  used  to  order  the  spectrum,  the  e.v.  are  C,,.  The  system  must  be 

solved  for  Cn  and  the  eigenfunctions  un,  vn,  wn,  pn  given  V,  c/a.  Re,  m  and  k; 

t  enters  only  through  V.  If  the  system  were  self-adjoint  the  Index  n  would  be 
the  radial  mode  numoer.  For  this  non-self-ad joint  system  it  is  unclear  how  to 
define  a  mode  In  a  general  and  unambiguous  manner.  A  mode  can  be  Identified 
for  large  t  by  calculating  the  e.v.  for  t  ♦  »,  where  the  radial  mode  number  is 
known  for  the  solid  rotation  solution,  and  then  tracked  as  t  decreases.  This 
Identification  Is  unambiguous  If  the  critical  layer  does  not  exist. 


8.  The  Numerical  Method. 


For  V  (r;t)  obtained  from  the  finite  difference  solution  to  (2.1),  the 
system  (4.2),  (4.4)  and  (4.5)  must  be  integrated  for  0  <  r  <  1.  The  numerical 

method  used  to  solve  the  e.v.  problem  is  a  shooting  method  with  iteration. 
The  numerical  integration  cannot  start  at  r  =  0  because  (4.2)  have  a  singular¬ 
ity  there.  Three  regular  linearly  independent  solutions  were  obtained  that 
satisfy  (4.5)  for  m  =  1  using  power  series  expansions  near  r  =  0  and  express¬ 
ing  V  as  a  power  series  in  r.  These  solutions,  evaluated  at  r  =  e,  are  used 
as  initial  conditions  for  the  numerical  integration  of  (4.2)  for  e  <  r  <  1; 

typically  e  *  .001.  Further  details  are  given  in  Reference  2. 

Integration  of  (4.2)  is  not  straightforward  because  the  coefficient  of 

the  highest  order  derivative,  Re~^,  is  small  for  the  values  of  Re  of 
interest.  To  insure  that  the  linearly  independent  solutions  remain  so  as  the 
integration  proceeds,  orthonormal ization  was  applied;  a  modification  of 
Davey's  technique9  was  used;  see  Reference  10  for  details.  Typically, 
orthonormalization  was  applied  at  50  equally-spaced  points  with  an  integration 
interval  of  .001;  these  parameters  are  varied  depending  on  the  case 
considered. 

Before  the  integration  process  is  started  a  value  of  C  must  be  specified, 
the  first  guess  for  Cn.  A  linear  combination  of  the  3  linearly  independent 

solutions  is  tested  to  see  if  (4.4)  is  satisfied.  This  test  requires  that  the 
characteristic  determinent  Z(C)  =0.  If  Z  *  0  within  a  certain  tolerance,  a 
new  value  of  C  is  obtained  using  Newton's  method;  Muller's  method  has  also 
been  used.  Equation  (4.2)  is  integrated  again  with  the  new  value  of  C,  etc. 
The  iteration  continues  until  a  convergence  test  on  the  iterates  is  satis¬ 
fied.  The  e.v.  is  then  known  and  the  constants  in  the  linear  combination  are 
determined  to  give  the  eigenfunctions.  One  of  these  constants  is  set  equal  to 
unity  for  convenience,  which  normalizes  the  eigenfunctions. 

Usually  a  spin-up  e.v.  history  is  required.  The  computation  is  started 
at  large  t  where  the  first  guess  can  be  obtained  from  the  Cn  for  solid  body 

rotation  or  the  inviscid  approximation  to  it.  For  smaller  t  the  first  guess 
is  obtained  by  extrapolation.  As  t  decreases  this  first  guess  may  not  be 
sufficiently  close  to  the  desired  e.v.  and  the  iteration  process  will  either 
converge  to  some  other  e.v.  or  diverge;  a  more  detailed  searching  process  is 
then  required.  In  our  computations  of  Cn  this  has  occurred  for  n  *  2  and 


9.  A.  Davey ,  " A  Simple  Numerical  Method  for  Solving  Orr-Soimerfeld  Problems," 
Quarterly  Journal  of  Mathematics  and  Applied  Mechanics ,  Vol.  26,  Part  4, 
1973 ,  pp.  401-411. 

10.  C.  U.  Kitchens,  Jr.,  N.  Gerber,  and  R.  Sedney ,  " Oscillations  of  a  Liquid 
in  a  Rotating  Cylinder:  Part  I.  Solid-Body  Rotation ,"  US  Amy  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  MD .,  BRL  Technical  Report 
ARBRL-TR-02081,  June  1978.  (AD  A0677S9) 


small  t  in  which  case  the  solution  of  the  e.v.  problem  becomes  more  labori¬ 
ous.  The  need  for  a  more  detailed  search  is  indicated  by  the  shape  of  the 

surfaces  Z(C)  and  will  be  illustrated  by  an  example. 

For  Case  1,  sections  of  the  surfaces  Z(C)  will  be  shown  for  C  in  the 
neighborhood  of  C£,  at  t  =  155.  For  illustrative  purposes  it  is  convenient  to 

plot  the  single  surface  |z|(CR,  Cj)  but  only  sections  of  that  surface  with 

planes  Cj  =  constant  are  shown  in  Figure  2.  For  the  parameters  of  Case  1,  m  = 

1,  k  =  3  and  t  =  155,  CR2  =  *072614  and  C^  =  .034266.  The  shape  of  the 

surface  changes  rapidly  in  the  neighborhood  of  that  point  for  small  changes  in 
CR  and  Cj.  The  first  guess  must  be  in  an  interval  aCr  *  .0005  and  aCj  = 

.0003  about  the  e.v.,  that  is  within  about  1%  of  C2I  otherwise  the  iteration 

converged  to  Ci.  The  high  peak  near  the  zero  of  |Z|  is  symptomatic  of  this 

behavior.  A  method  for  handling  this  situation  was  developed.  Note  that,  for 
n  =  1,  but  all  other  parameters  the  same  as  above,  the  |Z|  vs  Cp  curves  are 

monotonic  on  either  side  of  the  e.v.  and  the  first  guess  can  be  chosen  in  a 
much  larger  interval;  in  fact  extrapolation  is  sufficient. 

The  time  to  calculate  one  C  depends  on  many  parameters.  For  straight¬ 
forward  cases,  typical  CPU  times  for  one  C  calculation  are  1  minute  on  the  CDC 
7600  and  5  minutes  on  the  VAX. 


V.  THE  CRITICAL  LAYER 

The  Orr-Sommerfeld  equation  that  governs  the  perturbations  on  a  2-D  shear 
flow  can  be  regarded  as  a  prototype  for  the  system  (4.2).  The  former  is  4th- 
order  rather  than  6th-order  and  there  is  no  singularity  at  the  origin  in  the 
0-S  case.  However,  the  stiff  nature  of  the  0-S  equation,  because  of  the 

Re-1  factor,  and  the  possible  existence  of  a  critical  layer  make  it  a  useful 
guide  in  the  present  problem.  A  recent  review  of  various  aspects  of  the  criti¬ 
cal  layer  in  the  0-S  case  was  given  by  Stewartson. 11  A  discussion  of  it  using 
the  classical  solutions  to  the  0-S  equation  can  be  found  in  Schlichting1* 
together  with  a  description  of  the  famous  Schubauer  and  Skramstad  experiment 
which  verified  the  change  in  phase  of  the  streamwise  component  of  velocity 
across  the  critical  layer.  The  change  in  phase  of  the  Reynolds  stress  across 
the  critical  layer  is  essential  to  an  understanding  of  the  instability  of  a 
shear  flow. 

The  definition  of  the  critical  layer  is  best  appreciated  by  considering 
the  inviscid  perturbation  equations.  In  the  limit  Re  +  •  the  0-S  equation 
reduces  to  the  Rayleigh  equations  for  inviscid  disturbances.  For  a  neutral 


11.  K.  Stewartaon ,  "Marginally  Stable  Inviecid  Flows  with  Critical  Layers  " 
Journal  of  Applied  Mathematical  Vol.  27 t  1981 ,  pp.  133-175. 

12.  H.  Schliohting ,  Boundary  Layer  Theory ,  4th  Edition ,  McGraw-Hill  Book 
Company ,  New  Yorkt  1960,  Chapter  XVI. 
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disturbance,  if  there  is  a  point  at  which  the  phase  velocity  is  equal  to  the 
basic  flow  shear  velocity  (and  the  curvature  there  is  not  zero),  this  point  is 
called  a  critical  level  and  its  neighborhood  the  critical  layer.  The  critical 
level  is  a  singular  point  of  the  Rayleigh  equation;  a  discontinuity  in  the 
perturbation  velocity  is  implied,  violating  the  small  disturbance  assumption. 
There  are  ways  of  circumventing  the  singular  behavior.11  For  the  0-S  equation 
the  critical  level  is  a  turning  point  but  not  a  singular  point.  Some  of  the 
effects  of  the  critical  layer  are  different  for  amplifying  and  decaying 
disturbances.  The  latter  are  analogous  to  the  effects  in  the  present  work. 

The  inviscid  limit  of  the  perturbation  equations  for  the  rotating  fluid 

case  is  obtained  by  setting  Re-1  =  0  in  (4.2).  The  order  of  the  system  is 
reduced  from  six  to  two  and  this  system  is  analogous  to  the  Rayleigh  equa¬ 
tions.  Consider  a  neutral  disturbance,  i.e.,  Cj  =  0.  The  coefficient  of  the 

highest  order  derivative  contains  M  =  C  -  mV/r  and  the  C  is  real.  If  M  =  0 
has  a  real  root,  rc,  where  0  <  r  <  1,  the  equation  has  a  singular  point  at 

the  critical  level  rc  and  the  small  disturbance  assumption  is  violated.  The 

neighborhood  of  rc  is  called  the  critical  layer.  The  physical  interpretation 
of 

CK  =  mV/r  (5.1) 


at  r  =  rc  is  that  the  wave  frequency  is  an  integral  multiple  of  the  angular 

frequency  of  the  basic-flow,  indicating  a  resonance.  If  m  =  0  or  if  m  and 
Cr  have  opposite  signs  there  is  no  critical  layer.  The  nature  of  the  V  vs  r 
curves  shown  in  Figure  1  shows  that  an  rc  always  exists  for  small  t  if  m  *  0 

and  sgn  m  =  sgn  CR  but  will  not  exist  for  large  t.  There  can  be  a  critical 
layer  for  each  n. 

For  viscous  perturbations  (4.2)  do  not  have  a  singularity  when  M  =  0. 
The  r  for  which  M  =  0,  in  general  complex,  is  a  turning  point  of  (4.2). 
However,  the  neighborhood  of  real  rc  obtained  from  (5.1)  is  still  called  the 

critical  layer.  One  practical  consequence  of  the  existence  of  the  critical 
layer  is  that  the  eigenfunctions  can  develop  high  frequency  oscillations  of 
large  relative  amplitude.  These  occur,  most  notably,  for  small  t  and  large 
Re.  The  integration  scheme  and  the  number  of  significant  figures  in  the 
computation  must  be  capable  of  resolving  these  in  order  to  get  a  solution  to 
the  e.v.  problem  for  (4.2). 

Lynn13  showed  that  the  thickness  of  the  critical  layer  is  0(Re~^)  as  it 
is  in  the  0-S  case.  The  critical  layer  and  the  boundary  layer,  0(Re-1^),  can 
merge  at  the  lower  Re.  Stewartson11  showed,  analytically,  for  the  0-S  case. 


13.  Y.M.  Lynn ,  "Free  Oscillations  of  a  Liquid  During  Spin-Up /'  US  Army 
Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground ,  Maryland.  ARBRL 
Report  No.  1663 ,  August  1973.  (AD  A769710) 


that  the  interval  over  which  the  eigenfunctions  are  "violently  oscillatory"  is 
proportional  to  Cj  for  a  given  basic  shear  flow;  this  could  not  be  verified  by 

the  results  for  the  rotating  fluid  case. 


VI.  RESULTS 

Some  results  for  eigenfunctions,  e.v.,  and  the  critical  layer  will  be 
shown.  It  is  convenient  to  designate  the  three  wave  numbers  by  the  triplet 
(k,  n,  m)  corresponding  to  the  wave  numbers  for  the  (z,  r,  e)  directions. 


Table  1.  Effect  of  Neglecting  Diffusion. 


WITH 

1  DIFFUSION 

WITHOUT  DIFFUSION 

t 

CR 

Cj  x  103 

rc 

CR 

Cr  x  103 

rc 

rf 

600 

.149 

17.1 

.423 

.232 

27.2 

.668 

.618 

1245 

.0855 

3.14 

— 

.132 

4.44 

.392 

.368 

In  Figure  3  the  eigenfunction  Real  (w)  =  wR  vs  r  is  shown  at  three  times 
for  Re  =  5  x  105,  A  =  2.679,  ts  =  3,788  and  mode  (5,1,1);  laminar  Ekman  layers 
are  assumed.  In  Figure  3a,  t  =  7,000,  t/ts  =  1.85,  CR  =  8.354  x  10”2,  Cj  = 
8.363  x  10"*  and  there  is  no  critical  layer;  the  variation  of  wR  through  the 
boundary  layer  can  be  barely  discerned  on  this  scale.  In  Figure  3b,  t  =  1000, 
t/t$  =  .26,  CR  =  7.060  x  10”2,  and  rc  =  0.44  as  indicated  by  the  arrow  in  the 
figure;  the  rapid  variation  of  wR  in  the  neighborhood  of  rc  is  typical  of  the 
effects  of  the  critical  layer  on  the  eigenfunctions  for  Re  >  10*,  approxi¬ 
mately.  In  Figure  3c,  t  *  400,  t/ts  =  .11,  CR  =  3.015  x  10"2,  Cj  =  8.066  x 
10-2,  and  rc  =  0.66;  outside  the  critical  layer  wR  is  not  zero  although  it 
appears  to  be  on  this  scale.  The  oscillations  of  wR  are  centered  at  r  ■  rc. 
The  max  jwR|  is  greater  in  Figure  3c  compared  to  that  of  Figure  3a  by  a  factor 
of  1900.  For  the  conditions  of  Figure  3c,  if  a  =  628  rad/sec  (100  Hz),  £  = 
0.64  sec. 


Time  histories  for  Cj  and  are  shown  in  Figure  4  for  Case  1  and  modes 

(3,1,1)  and  (3,2,1).  The  CR1  curve  has  a  shallow  maximum  and  minimum  for  180 

<  t  <  220  but  on  the  scale  of  this  figure  it  appears  to  be  constant.  At  t  » 

45  it  has  a  maximum  which  is  typical  for  all  CR1  vs  t  curves;  note  that  CR  =  0 

for  t  =  0.  The  signficance  of  a  maximum  in  the  curve  is  related  to  a 

necessary  condition  for  projectile  instability.  If  the  nutational  frequency 
is  less  than  the  maximum  of  CR,  there  are  two  times  at  which  instability  might 

develop.  The  Cjj  and  CR2  curves  are  rather  typical,  but  the  Cj2  curve  is  not 

because  it  has  a  maximum  at  t  =  195  and  an  inflection  point  at  t  =  170;  these 
make  it  difficult  to  obtain  a  first  guess  for  C2.  The  problems  encountered  in 

calculating  C2  for  t  <  200  were  discussed  in  Section  IV-2.  Since  Cj2  <  Cjj 

for  t<140  the  n  =  2  mode  could  be  more  significant  than  the  n  =  1  mode  for 
projectile  instability. 

The  variation  of  rc  with  t  is  shown  in  Figure  5  for  both  Cases  1  and  2 

and  n  =  1  and  2.  Note  the  two  time  scales.  For  both  cases  rc*0  at  an  earlier 

time  for  n  =  1.  The  rc's  for  n  =  1  and  2  are  equal  when  CR1  *  CR2  but  Cjj  * 
Cj2  and  the  eigenfunctions  are  distinct.  For  Case  1,  n  =  1,  rc  =  0  at  t  =  225 
=  .489  ts  and  for  Case  2,  n  ■  1,  rc  »  0  at  t  «  9950  =  3.76  tst  which  gives  t  = 

13.2  sec.  Thus  for  Case  2  the  critical  layer  exists  over  a  substantial  part 
of  the  projectile  flight  time;  however,  its  effects  are  not  great  when  rc  is 
small. 

As  a  digression  from  the  presentation  of  results,  consider  the  implica¬ 
tions  of  using  (2.3),  which  neglects  diffusion  terms,  to  determine  the  basic 
flow  V  profiles.  Its  solution  is  (2.6)  for  r  >  rf  and  zero  otherwise.  In 

Table  1  the  e.v.  and  rc  determined  this  way  are  compared  with  those  using  the 
V  including  diffusion  for  Re  =  39,771,  A  =  3.12,  t$  =  1,245,  mode  (3,1,1) 

and  k  =  0.5.  Neglecting  diffusion  gives  large  errors  in  C  and  rc  and  the  rc 

is  never  zero  except  in  the  limit  t  +  »  . 

The  CR  time  history  for  Case  2,  mode  (5,1,1)  using  the  turbulent  Ekman 

layer  compatibility  condition  (2.5)  is  given  in  Figure  6;  the  maximum  occurs 
for  t  <  1,200.  Calculations  for  this  case  and  some  related  ones  were  used  to 
plan  projectile  firings  and  then  to  analyze  the  results  in  Reference  14.  In 


14.  W.P.  D'Amico ,  Jr. ,  "Flight  Data  on  Liquid-Filled  Shell  for  Spin-Up 
Instabilities ”  US  Ballistic  Research  Laboratory ,  Aberdeen  Proving 
Ground ,  Maryland ,  BRL  Memorandum  Report  ARBRL-MR-03334 ,  February  1984. 
(See  also  AIAA  Paper  83-2143 ,  August  1983.) 


planning  such  tests  the  sensitivity  of  Cp  to  variations  of  parameters  must  be 

considered.  The  results  for  two  such  variations  are  shown  in  Figure  6:  Re  was 
decreased  by  a  factor  of  10  keeping  A  fixed  which  decreased  Cp  and  A  was 

increased  by  10%  keeping  Re  fixed  which  increased  Cp;  the  relative  changes 

vary  with  t.  In  the  first  variation  Re  is  in  the  transition  range;  turbulent 
Ekman  layers  were  assumed.  There  is  a  large  change  in  rc(t)  for  the  first 

variation,  e.g.,  rc  =  0  at  t  =  3200  compared  to  9950  for  Case  2.  For  the 

second  variation  rc(t)  does  not  differ  much  from  Case  2,  e.g.,  rc  =  0  at  t  = 

11,700.  The  spin-up  times  are  1671  and  2912  for  the  first  and  second  varia¬ 
tions,  respectively. 

The  rapid  change  in  phase  of  the  velocity  across  the  critical  layer  in 
the  0-S  case  was  mentioned  earlier.  In  the  rotating  fluid  case  the  change  in 
phase  of  u'  for  0  <  r  <  1  will  be  shown.  From  (4.1) 

-C  t 

u'  =  (uR2  +  u,2)1'2  e  1  x  cos  Kz'  sin  [3  -  (CRt  -  mo)]  (6.1) 

tan  3  =  -uR/Uj  . 


The  phase  angle  3(r)  is  plotted  in  Figure  7  for  Case  2.  The  rapid  change  in 
phase  in  the  critical  layer  is  evident  for  t  =  2000,  6000,  and  9000;  the  mag¬ 
nitude  of  the  phase  change  in  the  critical  layer  is  approximately  223°,  180°, 
and  180°,  respectively.  Another  rapid  change  in  B  takes  place  in  the  sidewall 
boundary  layer.  This  phase  angle  change  is  a  sensitive  indication  of  the 
critical  layer.  As  t  increases,  the  r  interval  over  which  3  is  essentially 
constant  increases.  For  the  conditions  of  Figure  7,  3  =  constant,  except  in 
the  boundary  layer,  for  t  >  10,000  for  which  rc  =  0. 

The  wave  surface  u‘  =  constant  can  be  obtained  from  (6.1).  For  fixed  z* 
and  t  this  surface  appears  as  a  curve  in  the  (r,e)  plane.  For  certain  values 

2  2  1/2 

of  the  parameters,  the  variation  with  r  of  (uR  +  Uj  )  can  be  neglected  and 
the  wave  surface  is  given,  approximately,  by 


0  =  constant  -  3(c) 

This  curve  is  plotted  in  Figure  8  for  Case  2,  n  -  1  at  t  =  6,000.  If  there 
were  no  boundary  layer  or  critical  layer,  the  curve  would  be  a  straight  line 
through  the  origin. 

For  n  =  2  the  3(r)  curves  include  a  phase  change  across  the  critical 
layer,  as  for  n  =  1,  and  that  arising  from  the  fact  that  up  and  uj  must  each 

have  a  zero  since  the  radial  mode  number  has  increased  by  one.  In  Figure  9 
3(r)  is  shown  for  Case  2,  n  =  2  at  t  =  6,000  for  which  rc  =  .303.  The  total 

change  in  phase  for  0  <  r  <  1  is  about  400°,  considerably  larger  than  that  for 
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n  =  1  shown  in  Figure  7.  (That  this  is  not  always  the  case  is  shown  by  the 
next  example.)  In  addition  to  those  in  the  critical  layer,  uR  and  uj  have 

zeros  at  r  =  .420  and  .503,  respectively.  For  no  critical  layer,  say  t  * 
16,000,  there  is  an  abrupt  change  in  phase  of  180°,  rather  than  a  constant  6 
as  for  n  =  1. 

The  slopes  of  the  s(r)  curves  in  the  critical  layer  are  not  as  large  for 
smaller  Re.  In  Figure  10  0(r)  is  plotted  for  Case  1  at  t  =  160  for  n  =  1,  2 
for  which  rc  =  .304,  .176,  respectively.  Here  the  change  in  phase  for  n  =  2 
is  less  than  that  for  n  =  1. 


VII.  DISCUSSION 

The  theory  and  a  method  for  the  solution  of  the  spin-up  e.v.  problem  were 
presented  here.  The  theory  for  the  perturbed  flow  is  a  linear  one  but  not  for 
the  basic  flow.  Because  viscous  effects  are  important  in  the  boundary  layer 
at  the  cylinder  wall  and  in  the  critical  layer,  viscous  perturbations  are 
used.  The  theory  has  limitations  because  of  the  various  assumptions  that  are 
made.  The  theory  and  method  are  successful  in  the  sense  that  they  provide 
results  that  are  physically  meaningful  and  do  not  violate  intuition  or  the 
“physics  of  the  problem."  Other  investigators  have  worked  on  this  problem 
without  success  in  that  sense.  The  more  important  gauge  of  success  is  vali¬ 
dation  by  comparison  with  either  experimental  results  or  numerical  simulation. 

For  the  solid  rotation  case  and  m  =  0,  a  detailed  validation  was  given  in 
Reference  3  using  numerical  simulation.  This  provides  confidence  in  the 
treatment  of  (4.2),  and  the  solid  rotation  endwall  correction,  for  that 
case.  Spin-up  and  m  =0  could  be  validated  in  the  same  way;  it  would  require, 
relatively,  more  analysis  to  reduce  the  numerical  data.  At  the  present  time 
there  is  no  numerical  simulation  to  validate  the  m  =  1  case. 

To  the  authors'  knowledge,  the  only  reported  measurements  of  C  for  m  =  1 
are  in  Reference  6.  Some  experimental  results  for  Cr  vs  t,  using  the  methods 

of  Reference  6,  are  shown  in  Figure  11;  Stergiopoulos,  private  communication. 
The  fact  that  the  data  are  in  groups  of  3  points,  is  an  artifact  of  the  data 
reduction  process.  For  each  group  Re  =  constant;  for  all  groups  4.39  <  Re  x 
-4 

10  <  4.45,  a  negligible  variation.  The  aspect  ratio  A  =  0.600  and  the  mode 

is  (1,2,1).  The  scatter  in  the  data  is  about  ±3%  except  at  t=245  where  it  is 
±5%.  The  calculated  Cr,  using  <  =  0.443,  is  within  the  scatter  of  data.  This 

comparison  validates  the  calculation  of  Cr  for  these  parameters  over  the 

range  of  the  data  0.724  <  t/t$  <  2.252.  For  large  t,  essentially  at  solid 

rotation,  the  experimental  and  calculated  results  differ  by  1.5%. 6  The  rc(t) 

curve  shows  that  the  critical  layer  exists  over  a  considerable  range  for  which 
data  is  presented. 

The  percent  differences  of  CR,  the  frequency  with  respect  to  an  inertial 

frame,  are  considerably  larger  than  those  for  the  frequency  with  respect  to 
the  rotating  frame,  Cr_.  For  solid  rotation,  or  steady  state,  it  is  more 


conventional  to  work  in  the  rotating  frame.1*  The  relationships  between  the 
two  frequencies  and  their  percentage  changes  are 


CR  =  1  -  cRr 


ACr/Cr  =  -UCRr/cRr)  (1  -  CR)/CR 

Comparisons  of  experimental  and  calculated  CR  for  modes  (2,1,0)  and 

(1,1,1,)  give  the  same  conclusions  as  for  the  (1,2,1).  Comparisons  of  Cj  are 

not  shown  for  two  reasons:  (i)  the  experimental  error  in  the  determination  of 
Cj  can  be  quite  large  and  (ii)  not  accounting  for  the  Ekman  layers  in  the 

calculation  gives  Cr's  about  one-half  of  the  proper  values;  the  effect  on  Cn 
is  1  -  2%. 

There  are  no  data  in  the  spin-up  range  from  projectile  firings  that  can 
be  used  to  validate  this  theory.  Some  indirect,  qualitative  results  are  given 
in  Reference  14  which  show  consistency  with  the  theory  presented  here. 
Further  validation  will  probably  depend  on  laboratory  experiments. 


ACKNOWLEDGMENT 

The  authors  wish  to  thank  Ms.  Joan  M.  Bartos  for  her  work  on  a  large 
number  of  tasks  which  made  the  programs  used  here  operational.  Her  pro¬ 
gramming  and  editing  were  essential  to  the  completion  of  this  work. 


22 


- Re- 4.97*  103,  A -3.30,  t,  -  460 

- Re  =  1.99  *  10‘.  A  -  5.20,  ♦„  -  2648 


t-  1350 
t=  4050 
1-  13500 


//; 

//h 


yy/A'I 


/  /  / 

/A  / 


» =  39.9 
1=  131.2 
I  =  262.5 


Figure  1.  V  vs  r  for  Cases  1  and  2  at  Three  Times;  k  = 
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Figure  2.  The  Magnitude  of  the  Characteristic  Determinant,  |Z|,  vs  Cr 
Three  Values  of  Cj  for  Case  1,  m  *  1,  k  *  3,  t  »  155.  Cr?  * 
.072614,  Ci2  *  .034266. 


Figure  4.  Time  Histories  for  Cj  and  Cg  for  Case  1,  Modes 
(3,1,1)  and  (3,2,1),  ts  =  460,  k  =  0.5. 


so 


R««  1.989  *  10*.  A  -  5.20 
k  -  5,  n  -  1 


Figure  7.  The  Phase  Angle  of  u‘  for  Case  2,  Mode 
(5,1,1)  at  Three  Times. 


Figure  8.  In  the  (r,e)  Plane,  the  Approximate  Wave  Surface  u'  *  Constant  for 
Fixed  z'  and  t  for  Case  2,  Mode  (5,1,1)  at  t  =  6,000.  The  Circular 
Arc  Portion  Is  at  the  Critical  Level  r  =  r  =  .20;  the  Change  in  e 

Near  r  *  1  Occurs  In  the  Boundary  Layer. 


Figure  9.  The  Phase  Angle  of  u'  for  Case  2,  Mode  (5,2,1)  at  t  = 
The  Critical  Level  is  at  rc  =  .303. 


r 


Figure  10.  The  Phase  Angle  of  u'  for  Case  1,  Modes  (3,1,1)  and 
(3,2,1)  at  t  *  160. 


Figure  11.  CR  vs  t:  Experimental  Points  and  Calculated  Values 


r  vs  t,  k  =  0.443. 
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